Impurity Doping Effect in High T c Superconductors 



Yong-Jihn Kim and K. J. Chang 
Department of Physics, Korea Advanced Institute of Science and Technology, 

Taejon 305-701, Korea 



Abstract 



It has been observed that impurity doping and/or ion-beam- induced dam- 
age in high T c superconductors cause a metal-insulator transition and thereby 
suppress the critical temperature. Based on our recent theory of the weak lo- 
calization effect on superconductivity, we examine the variation of T c with 
increasing of impurity concentration (x) in Lai.85Sro.i5Cui_ x A x 04 systems, 
where A = Fe, Co, Ni, Zn, or Ga. We find that the doping impurity de- 
creases the scattering matrix elements for electron-electron attractions, such 
as V nn / = —V[l — v ^ (} n{L / '£)], where L and I are the inelastic and elastic 
mean free paths, respectively. Using the mean free path I determined from 
resistivity data, we find good agreements between our calculated values for 
T c and experimental data except for Ni-doped samples, where Ni impurities 
may enhance the pairing interaction. 



PACS numbers: 74.20.-z, 74.40.+k, 74.60.Mj 
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Although impurity scattering may provide a clue to understanding the superconducting 
mechanism of high T c superconductors, its effect is not well understood yet. 1-3 In strongly 
correlated systems such as high T c superconductors, impurity potentials are significantly 
renormalized by strong electron-electron interactions. 4 In particular, (unitarity) resonant 
impurity scattering seems to be important. Using the marginal Fermi liquid theory, Kotlier 
and Varma 4 showed that scattering by ordinary impurities is close to the unitarity limit. 
Recently, it was argued that very low concentrations of impurities in non-Fermi-liquid sys- 
tems lead to a vanishing of density of states at the chemical potential and infinite resistivity 
for T — > 0. 5 In a one-channel Luttinger liquid with repulsive interactions, Kane and Fisher 6 
showed that electrons are completely reflected even by the smallest scatterer as if it were 
effectively infinite. Hirschfeld and Goldenfeld 7 considered the effect of resonant scattering 
on the penetration depth of a rf-wave superconductor. Poilblanc, Scalapino and Hanke 8 
showed by an exact diagonalization technique for small clusters that strong correlations give 
rise to resonant scattering due to local defects in the t — J model. In addition, Ziegler and 
his coworkers 9 introduced a simple static approximation to the effective impurity potential 
which reproduces resonant scattering for both two-dimensional t — J and Hubbard models. 
Nevertheless, the impuirty effect on the transition tempearture of high T c superconductors 
remains to be unsettled, because the effect of resonant scattering on the superconducting 
temperature is not much different from that of Born scattering. 3 ' 10 

We note that experimental results show consistently that impurity doping and ion-beam- 
induced damage give rise to a metal-insulator transition and thereby suppress T c n - 15 Re- 
cently, we studied the weak localization effect on low T c superconductors, 16 and found that 
for two-dimensional systems weak localization decreases effective interactions, similarly to 
conductivity. Here we use the same theoretical approach to investigate the impurity doping 
effect on the superconducting temperature for Lai. 85 Sr a i 5 Cu04 with Cu replaced with 3d 
transition elements or Ga. The effect of strong correlation is taken into account by the 
simple renormalization of the impurity potential. 9 The mean free path I is determined from 
resistivity data. We find good agreements between our theoretical calculations and experi- 
mental data except for Ni-doped samples. In Ni-doped case, the suppression of T c is much 
slower than the theory predicts, implying that Ni may enhance the pairing interaction. 

Let us consider a simple pairing potential, V(|ri — r 2 |), between electrons of a pair. The 
matrix elements of V(r) for the plane wave states are 

V kk> = J j dr l dr 2 e^ k '^-^V{\r l -v 2 \)e it ^- r ^ 

= v + v x k ■ y + ^[3(fc • Vf - 1] + • • • . (i) 

For a s-wave pairing, V^p is assumed to be a constant Vb < 0, whereas for a ci-wave pairing, 
Vj££, depends on the unit vectors k and k , such as 

%' = y[3("') 2 -i]- (2) 

In the presence of impurities, the pairing interaction between scattered basis pairs (ip n ,ipn) 
and (ip n ',ipn') is given by 17 

V nn , = JI dr 1 dr 2 ^(r 1 )^(r 2 )^(|r 1 - r 2 |)^(r 2 )^(ri), (3) 
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where denotes the time-reversed partner of ip n . The metal-insulator transition driven by 
the impurity doping may be understood by the localization of wave function. 18 ' 19 In high T c 
superconductors, resistivity data showed that correlation enhances significantly the strength 
of impurity potential. Kaveh and Mott 19 derived the weakly localized scattered states in the 
form of power-law and extended wave functions for two-dimensional systems, 

Mr) = A 2 e rk - r + B 2 — , (4) 



Aj = l-2nB 2 2 ln(L/t), B\ = -^— £ , (5) 

where £ and L are the elastic and inelastic mean free paths, respectively. We can calculate 
V nn i by substituting Eq. (4) into Eq. (3). Since the Cooper pairs in a BCS condensate 
form bound states, only the power-law wave functions within the BCS coherence length are 
relevant. 20 Thus, the Cooper pair wave functions basically consist of the plane waves with 
reduced amplitudes. The resulting matrix elements are 21 ' 22 



(6) 



for both the s- and rf-wave pairings. In the ci-wave pairing, Eq. (6) does not include the 
effect of impurity scattering in the Borm limit. Then, the effective coupling constant X e ff 
becomes 

\ eff = \[l--^- £ ln(L/e)}, (7) 

which satisfies the modified BCS T c equation, 

T c = 1.13e c exp(-l/A e// ). (8) 

Xiao and his coworkers 11 conducted a systematic study on the effect of Cu-site doping 
in Lai.85Sro.i5Cu04 polycrystals and determined the relation between T c and dopant level 
for Fe-, Co-, Ni-, Zn-, Ga-, and Al-impurities. Cieplak and his coworkers 12 also measured 
resistances for the same materials with different impurities (Fe, Co, Ni, Zn, and Ga) down 
to 50 mK and found the metal-insulator transitions produced by impurities. Here we use 
Eq. (8) to calculate T c as a function of dopant concentration, assuming that the cutoff 

energy e c = 450K and the Fermi wave vector Uf = 0.2A . In this case, A is chosen to give 
T co = 37K without impurities. The inelastic mean free path is obtained from disordered 
two-dimensional systems, which is employed for all the systems considered here, with the 
1/y/lnT dependence removed, i.e, L = SOOA/a/T- 18 Because of the Sr doping, the elastic 
mean free path is estimated to be about 50A without Cu-site doping. 13 Then the total elastic 
mean free path £ tot (in unit of A) is 

1 1 1 
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where £i mp is the mean free path caused by impurity doping for Cu. For two-dimensional 
layered systems such as high T c superconductors, the mean free path £ imp can be obtained 
from the Drude formula, kpt = h/R n e 2 , where R n is the sheet resistance, defined as R n = 
p/d. In this case, d is the interlayer distance (6.5 A in the LSCO system) 23 and p is 
resistivity From experimental data by Cieplak et al., 12 the resistivities of Ga-, Zn-, Fe-, 
Co-, and Ni-doped systems can be fitted such as pc a = 198.41 x c, pz n = 147.11 x c + 40, 
PFe = 227.4 x c + 33.7, pc = 190.67 x c, and p^i = 164.35 x c + 12, respectively, where c 
denotes the concentration of dopants. 

The superconducting transition temperatures are plotted as a function of dopant con- 
centration and compared with experiments in Figs. 1 and 2. For nonmagnetic (Ga and Zn) 
and magnetic (Fe and Co) impurities, the agreements between theory and experiment are 
satisfactory, considering the simple approximation of strong correlation effect. It is clearly 
seen that the impurity doping gives rise to a metal-insulator transition and suppresses T c . 
We find that only the strength of (effective) impurity potential is important and the effect of 
nonmagnetic impurities on destroying superconductivity is similar to magnetic impurities, 
which agrees well with experiments. Although the Zn and Ga atoms have filled d-shells 
whereas the Fe, Co, and Ni atoms have unfilled d-shells, it appears that the variation of T c 
with impurity doping does not depend significantly on the valence states of impurities. We 
find that the initial drops of T c for all the samples are larger than the theory predicts. 

Fig. 3 shows the variation of T c for Ni-doped samples. We find a large deviation between 
the theoretical and experimental curves. Although the origin of the discrepancy is not clear, 
it seems that Ni may enhance the pairing interaction. It is interesting to note that a Ni 
impurity in YBa2Cu307_ x has an unpaired spin of S — | rahter than S — 1, as expected 
for Ni 2+ . 24 It is also possible that Ni-doped samples may develop some inhomogeneous 
granular structures 25 or stripe phases 26 and their superconducting phases may persist for 
very high impurity concentrations. However, x-ray diffraction measurements in Ni-doped 
single-crystals YBa2Cu307_ x did not support this idea. 27 

In the underdoped regime, experiments showed that impurity scattering produces a large 
residual resistivity close to (or even larger than) the unitarity limit and the metal-insulator 
transition occurs almost simultaneously when T c drops to zero. 13,14 Consequently, our theory 
which considered the optimum doping regime is also applicable for the underdoped regime. 
As the hole density increases from the optimum doping, T c drops to zero more quickly before 
the metal-insulator transition is reached. This behavior may be related to the dimensional 
crossover from two-dimension to three-dimension. Nevertheless, our theory may explain 
the decrease of T c in moderately overdoped samples. Recently, Suryanarayanan and his 
coworkers 28 found the recovery of superconductivity in Y 1 _ x Ca x SrBaCu 2 .6Al .4O 6+z when 
Y is substituted by Ca. This result may be understood if we consider the importance 
of the metal-insulator transition caused by impurity doping. Since Ca increases the hole 
density in the Cu-0 planes, the mobility edge is shifted with respect to the Fermi energy, 
so that the system becomes metallic and superconducting. The transition between the 
insulating and superconducting phases was also explored in thin films such as DyBa 2 Cu 3 07, 29 
Nd2- x Ce x Cu04, 30 and YBa2Cu307, 31 where the two-dimensional localization is the main 
reason for the suppression of superconductivity and the critical sheet resistance lies in the 
range of 6 to 8 kQ. In YBa2Cu307__ x /PrBa 2 Cu307_ x superlattices, 32 ' 33 the decrease of T c 
in isolated YBa 2 Cu 3 07_ x layers with increasing sheet resistance is on the scale of h/e 2 , as 
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observed for ordinary superconductors. 32 This behavior may also be explained in terms of 
our theory. 

In conclusion, we have studied the impurity doping effect on the superconducting tem- 
peratures of Lai.85Sr . 15 Cui_ x A x O4 systems (A=Fe, Co, Ni, Zn, and Ga) based on our recent 
theory of weak localization effect on superconductors. The strong correlation effect is taken 
into account by a simple renormalization of the effective impurity potential. The calcu- 
lated values for T c are in good agreement with experiments, while in Ni-doped samples the 
agreement was poor, implying that Ni may enhance the pairing interaction. 
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Figure Captions 



Fig. 1 Variation of T c with dopant concentration for Lai. 85 Sr i 5 Cui_ x A x 04 

(A=Ga and Zn). Experimental data are from Ref. 11. 
Fig. 2 Variation of T c with magnetic dopant content for Lai.85Sr .i5Cui_ x A x O4 

(A=Fe and Co). Experimental data are from Ref. 11. 
Fig. 3 Variation of T c with Ni concentration for Lai.85Sro.i5Cui„ x Ni x 04. Experimental 

data are from Ref. 11. 
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Fig. 2 
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